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Abstract: We study relativistic hydrodynamics of normal fluids in two spatial dimensions. 
When the microscopic theory breaks parity, extra transport coefficients appear in the hy- 
drodynamic regime, including the Hall viscosity, and the anomalous Hall conductivity. In 
this work we classify all the transport coefficients in first order hydrodynamics. We then use 
properties of response functions and the positivity of entropy production to restrict the possi- 
ble coefficients in the constitutive relations. All the parity-breaking transport coefficients are 
dissipationless, and some of them are related to the thermodynamic response to an external 
magnetic field and to vorticity. In addition, we give a holographic example of a strongly 
interacting relativistic fluid where the parity-violating transport coefficients are computable. 
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1. Summary 

Hydrodynamics is an effective long-distance description of many classical or quantum many- 
body systems at non-zero temperature. The form of the hydrodynamic equations is dictated 
by the symmetries of the microscopic Hamiltonian, and is not sensitive to the precise nature 
of the short-distance degrees of freedom. When the microscopic description exhibits Lorentz 
invariance, the collective flow is described by the relativistic analogue of the Navier-Stokes 
equations. For normal fluids with an unbroken global U{1) symmetry (such as baryon num- 
ber), the hydrodynamic equations take the form [1,2], 

V^T'^'^ = F'^^ , V^J^ = 0. (1.1) 

Here T^'^ is the energy- momentum tensor of the fluid, is the U{1) symmetry current, and 
we have allowed for the possibility of coupling the fluid to an external non-dynamical gauge 
field (with field strength F^^^) and metric g^u (with covariant derivative V^). The gauge 
field couples to the conserved current J^. The relativistic analogue of the Navier-Stokes 
equations have a wide range of applications. For example, their study in 3-|- 1 dimensions has 
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led to significant progress in understanding the quark-gluon plasma [3]. In 2+1 dimensions the 
equations were proposed as an effective description of thermo-magnetic transport in cuprates 
[4] and in graphene [5]. 



A complete hydrodynamic description must, of course, supplement equations (1.1) with 
constitutive relations which express T^'^ and in terms of macroscopic parameters such as 
the local fluid velocity u^, local temperature T and local chemical potential ^. A conventional 
description of the constitutive relations can be found, for example, in the classic textbook by 
Landau and Lifshitz [1]. In this paper we take a closer look at the equations of relativistic 
hydrodynamics in 2 + 1 dimensions, and argue that the canonical constitutive relations need 
to be modified when the microscopic theory does not respect parity P. By parity we mean 
invariance under reflection of one of the spatial coordinates.^ An example of a P-violating 
system is the theory of interacting massive Dirac fermions in 2+1 dimensions; the mass term 
breaks parity. 

Our expressions for the constitutive relations may be written as follows: 

T^^'' = eou^'u" + (Po - CVau" - xbB - xn^) A^'^ - 7?a'^^ - r?^^'^ , (1.2a) 
= /)on^ + oV^" + + xeE^" + XT^^"'Pu^VpT . (1.2b) 

The tensor quantities appearing in the constitutive relations ( |1.2| ) are 

n = -e^^'Pu^V.Up, B = -^e^'-Pu^F.p, (1.3a) 

E^' = Pt'^u^, = E>' - TA'^'^V^^, (1.3b) 



and 



fl.3c) 



Ef" = e^'^Pu^Ep , = e^'^Pu^Vp , (1.3d) 

~^.u ^ 1 + e""''n,a/) . (1.3e) 

The thermodynamic parameters -Po(/^) T), eo(^, T) and po{n, T) are the values of the pressure, 
energy density and charge density respectively in an equilibrium configuration in which B = 
= 0, where B is the rest-frame magnetic field and Q the vorticity. They satisfy 

dPo = sodT + podfi , (1.4) 
€o = -Po + sqT + pon, (1.5) 

where sq is the entropy density. The velocity field is denoted u^^ and is normalized so that 
u^Ua = — 1. In this paper we study hydrodynamics to first order in derivatives. For counting 



^Parity can always be defined as reflection of one of the spatial coordinates. In a 3 + 1-dimensional theory, 
parity combined with a rotation is equivalent to a reflection of all three coordinates. In 2 + 1 dimensions, a 
reflection along both spatial coordinates is equivalent to a rotation. 
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purposes, the derivatives of and ^^t^ are of the same order as derivatives of the hydrody- 
namic variables. As a result, we take the magnetic field B and vorticity as first order in 
derivatives, and work to linear order in B and 0,. 



The remaining parameters in (1.2) characterize the transport properties of the fluid, or 
its thermodynamic response. The shear viscosity t/, bulk viscosity and charge conductivity 
a are the canonical dissipative transport coefficients and must satisfy 

rj^O, C^O, a^O, (1.6) 

as a consequence of either positivity of the divergence of the entropy current, or positivity of 
the spectral functions in the corresponding Kubo formulas. The Hall viscosity fj and a new 
parameter a, which contributes to the Hall effect in the absence of external magnetic fields, 
are both dissipationless. Our analysis does not constrain the values of a and fj,^ 

^eM, ctGR. (1.7) 

The remaining parameters xb, Xn, Xe and xt are not independent, and are specified in terms 
of three thermodynamic functions, M.b{T,ij-), M.Q{T,fj,) and /o(T), such that 

OMb p fdMn \ 
Xe = ^-Ro[-^-Mb), (1.8c) 

TxT = (T^ + -Mb)-Ro (T^ + + UT) -2Mn), (1.8d) 



where we have defined Rq = po/i^o + ^o)- All derivatives in ( |1.8| ) are evaluated at constant 
or T except for dPo/dcQ and dPo/dpo which are evaluated at constant po and eo respectively. 
The Kubo formulas for the parameters appearing in the constitutive relations (^]^) are 

7] = lim^-^{6ikSji - eifcej7)ImGj^''''(a;,0) , fj = hni ^^ifcEji Im G*^'^'(a;, 0) , (1.9a) 

a = lim -^dij Im G^f^{uj, 0) , cr + xe = lim -^eij Im G'a^{u, 0) , (1.9b) 

along with 

C = lim -3-(^,,-5fciImG)j''='(a;,0) , (1.9c) 



^When T = 0, it has been argued that f] should coincide with the angular momentum density of the ground 
state of the systems considered in [6,7]. 
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and 



eijk^ f dPp ^oo.oj ^ 



fc-s-o k^ V oer ^ 





(l.lOa) 




(1.10b) 


°(0,fc)), 


(1.10c) 


r°(o,^)), 


(l.lOd) 



XE = ^f^;-^[&R^^,k)-RoG'i 

TxT = I lim ^ (g^/°(0, fc) - RpG\ 

where ejj is an antisymmetric tensor with eyi = 1. Here Gr{uj, k) denotes the retarded Green's 
functions, 

in the thermal equilibrium state at = O and fl = 0, defined by varying the one-point 
functions with respect to the appropriate sources. One important difference between the 
Kubo formulas (1.9) and (I.IO) is that the former are given in terms of zero-momentum 
response functions, while the latter are given in terms of zero-frequency response functions. 
As emphasized in [8], the zero- frequency response functions are inherently Euclidean, and 
therefore only contain thermodynamic information. For this reason the parameters xb, Xn, 
Xe, XT are not transport coefficients, but should be thought of as thermodynamic quantities, 
consistent with {^^). We will refer to rj, fj, a and a as transport coefficients and to xb, Xn, 
Xe and xt as thermodynamic response parameters. 

Our parametrization of the constitutive relations ( |1.2D was not general in the sense that 
we have chosen a particular out-of-equilibrium definition of energy density, charge density, 
and fluid velocity. Such a choice is referred to as a "frame". The choice ( |1.2| ) is usually 
referred to as the Landau frame. We will find it convenient to use an alternative frame that 
is naturally suited to fluids whose thermodynamics depends on B and Q. Indeed, a static 
magnetic field B does not lead to an increase in (fluid) energy and therefore may be non-zero 
in equilibrium. Similarly, on a compact manifold, one may have non-zero vorticity and 
still remain in thermal equilibrium, e.g., a system which executes rigid rotation. For such 
equilibrium states the pressure is P = ^(T, /i, 0), so that 

dP dP 

dP = sdT + pdfi + —B + —n, (1.11) 
€ + P = sT + np. (1.12) 

Here and in the rest of this paper, all thermodynamic derivatives with respect to B and O are 
evaluated at B=0 and Q=0. The constitutive relations in a 'magnetovortical' frame which is 
adapted to the thermodynamic relation ( |1.11| ) are given by 

+ {P- CVan" -XBB-iun) AP" - rjaP" - ^a^'^ , (1.13a) 
J/^ = (p-MB^)u'' + aV'' + aV'' + XEE'' + XTeP''Pu^VpT, (1.13b) 
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where A4b and M.^ are 



dP dP 



The role of the undetermined function fn{T) is unclear. The expressions ( |1.8| ) and ( |1.14 ) 



then determine the parameters in the constitutive relations in terms of thermodynamic 

derivatives, 

dP dP , , 

de ^ fdP de \ dp ^ f dP dp\ 

We note that it is also possible to present our results in a frame-invariant form along the lines 
of the analysis carried out in [9]. We describe this in Section 0. 



The constitutive relations ( |1.2| ) (or (1.13)) together with the subsequent relations for the 
transport coefficients and thermodynamic response parameters are the main results of this 
paper. The relations ( |1.6| ) and the Kubo formulas for rj, a, and C are well known, while 
Kubo formulas for fj were discussed recently in [10]. Our strategy for obtaining the relations 
(|1.6D, (|1.7D, (|1.8[), and ( [1.10 ) involved the imposition of several physical constraints on the 



constitutive relations. These constraints amount to requiring that the response functions of 
a hydrodynamic theory must: (i) obey positivity constraints, (ii) have their zero-frequency 
limits coincide with the corresponding thermodynamic susceptibilities, and (iii) transform 
covariantly under time-reversal, T. In addition we ensure that (iv) a local version of the 
second law of thermodynamics holds. Some of these constraints are more familiar than 
others. Relations (i) and (iv) have often been used in the literature [1], and (iii) is the basis 
for Onsager's reciprocity relations [11,12]. 

Parity-violating systems in 2 + 1 dimensions have been considered in the condensed 
matter literature. The simplest such example is a theory of free massive Dirac fermions at 
zero temperature and in 2+1 dimensions. Parity breaking leads to a remarkable transport 
property at zero temperature: the Hall conductivity is quantized although no magnetic field 
is present [13]. This is an example of the anomalous Hall effect [14]. The transverse response 
to a thermal gradient (the thermal Hall conductivity) was recently discussed in several classes 
of topological insulators in 2+1 (and 3+1) dimensions, and related to anomalies in various 
dimensions [15]. In 2+1 dimensions the parity-odd analogue of the shear viscosity, which is 
called the Hall viscosity, has been studied from the condensed matter physics perspective in [6, 
16-19], using an effective field theory [7,20], and also using the AdS/CFT correspondence [10, 
21-23]. 

Parity- violating transport effects were also studied in 3+1 dimensions with QCD applica- 
tions in mind: using field theory techniques [24-28], in the context of hydrodynamics [29-32], 
and in parallel using the AdS/CFT correspondence along with hydrodynamics [33-41]. A 
relation linking parity-odd transport in 3+1 dimensions with the chiral anomaly was first 
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found in [42,43], and later effects of a gravitational anomaly were considered in [44,45]. Re- 
cently, parity-odd transport in superfluids was discussed in [9,46,47]. An effective field theory 
for non-dissipative transport in 1+1 dimensions was suggested in [48]. For considerations of 
hydrodynamics in arbitrary dimensions see [32,49]. 

The rest of this paper is organized as follows. In Section |2| we construct the most gen- 
eral constitutive relations allowed for a relativistic P-violating system in 2+1 dimensions. 
In Section ^ we construct an entropy current with positive divergence, and determine the 
ensuing constraints on the constitutive relations. We independently derive restrictions on 
these constitutive relations in Section ^ using linearized hydrodynamics, by computing the 
retarded Green's functions and imposing the conditions (i)-(iii) described above. We discuss 



an alternative hydrodynamic frame ( 1.13 ) in Section g, which provides a more transparent 
picture of 2 + 1-dimensional thermodynamics in the presence of non-zero B and fi. We check 
our results against an AdS / CFT computation in Section ^ and conclude our analysis with a 
discussion of the results in Section 0. 

2. The hydrodynamic expansion 

In hydrodynamics, the chemical potential, temperature and velocity field are allowed to vary 
slowly in space and time. The four equations of motion which determine the values of the 
hydrodynamic variables are energy-momentum conservation and charge conservation, while 
the explicit relations between the energy-momentum tensor (and the current) and the hydro- 
dynamic variables are called constitutive relations. Given a time-like vector (satisfying 
Ufj,u^ = — 1), the energy- momentum tensor and the current can be decomposed into pieces 
which are transverse and longitudinal with respect to u, 

T^u = £u^Uu + VA^,^ + {q^Uu+qyU^) + t^^, , (2.1a) 
J^,=Nu^ + j^, (2.1b) 

where, as before, A^j^ = g^u + u^Uy projects onto the space orthogonal to the velocity field. In 
this decomposition, "P, and M are Lorentz scalars, (?^, t^^^ and are transverse, u^g''^ = 0, 
u^t^^ = 0, u^j^ = 0, and t^^ is symmetric and traceless. Some readers may be familiar with 
a decomposition of the energy-momentum tensor and current into an ideal and dissipative 
part. For example, one may write 

T^^ = eQU^'u'' + PqA^^ + r'^^ , 

( 2.2 ) 



where Pq, eo and po where defined below (L3). We point out that the decomposition in ( |2.1D is 
of a different nature — it is a decomposition into scalar, transverse vector and transverse tensor 
modes which can be carried out for any tensor and vector. In other words, £, V and M do 
not necessarily take on their values in equilibrium. As a consequence, the scalars, transverse 
vectors and transverse tensors in ( |2.1aD and ( 2.1b| ) depend, a priori, on the hydrodynamic 
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variables fx, T and and on quantities built from their derivatives. Needless to say, one 
may easily go from ( |2.2| ) to (^]^) by comparing appropriate terms in the current or energy- 
momentum tensor. For example: M = po — u^T^. 

Out of equilibrium one can redefine the fields u^{x), T{x), and p{x) in a way that 
simplifies the decomposition This four-parameter field redefinition is referred to as a 

choice of frame [1] (for a detailed recent discussion see [50]). In what follows we will choose 
a conventional Landau frame, in which the four-parameter ambiguity is fixed by requiring 
that q'^ = and that £ and M retain their values in an equilibrium configuration with zero 
magnetic field and zero vorticity, i.e. £ = and M = po- This choice of frame gives 

J IX = PoUfi + J^l , 

where the transverse current j'^ and transverse traceless tensor tf^,^ vanish in the equilibrium 
state, and similarly 

V = Po + (.ttc«o°;J • (2.4) 

We now need to specify the constitutive relations which express the energy-momentum 
tensor and the current in terms of the hydrodynamic variables u^, T, and p, their derivatives 
and possible electromagnetic and gravitational sources. Since we consider small deviations 
from thermal equilibrium, we may expand V, t^y, and to first order in derivatives of 
hydrodynamic variables. As is generally the case in effective field theories, we must allow 
all possible one derivative contributions to V, t^^, and consistent with the symmetries of 
the system, but rule out those expressions that are forbidden by physical constraints such as 
thermodynamic laws, unitarity and time reversal symmetry. In the remainder of this section, 
we will classify all possible independent contributions to V, j^, and t^^'^ . By independent we 



mean expressions which are inequivalent under the equations of motion (LI) to first order in 
derivatives. The additional constraints which need to be implemented in order for the theory 
to satisfy all physical requirements will be described in Section ^ and Section ^ 

In formulating the first-order constitutive relations, we take the external fields to be 
small, with the field strengths F^y and the connection coefficients T^^u of the same order as 
gradients of the hydrodynamic variables. This is the scaling required to study the response of 
the fiuid to sources to first order in a derivative expansion. Magnetohydro dynamics or fluid 
dynamics with large values of vorticity would require a separate treatment. 

To carry out a classification of the scalars, vectors and tensors it is convenient to sup- 
plement the transverse projector A^^/, satisfying = A, with a transverse antisymmetric 
tensor S^jy = e^^ypU^ which satisfies Ti^^u'^ = 0, = —A and S • A = S. A generic vec- 
tor or pseudovector can be projected into orthogonal components in the plane transverse to 
via A^jy and Ti^u- This allows us to straightforwardly write down all possible structures 
contributing to the constitutive relations. 

At first order in derivatives there are three scalars, u^V^T, u^^V ^{p/T), I^^^'^V ^Uy = 
V^u^^, and two pseudoscalars, Tj'^'^V^Uy = —Q,, ^Ti^^Ff^y = —B that one could construct 
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out of the hydrodynamic variables. However, since there are two scalar equations of motion: 



and u^V^T^'^ = u^F^'^J^, only one of the three scalars is independent. We take 



this scalar to be A^'^V ^u^, = V ^u^. Thus, there is one scalar and two pseudoscalars which 



may contribute to V to first order in derivatives, 

V = Po-XBB-xn^ 



(2.5a) 



where Pq is the local thermodynamic pressure, and ^ is the bulk viscosity. The second and 



third terms in ( 2.5a ) are forbidden in parity- invariant systems, but are allowed once parity is 
broken. 

Next we consider the tensors. Since u^V^T and Uj^iW p{^/T) have no transverse projec- 
tions, it is sufficient to focus on projections of V fj,Uy. The two structures A^^^Vpu'^^ and 
T,P^f^\/ pu'^^ with circular brackets denoting a symmetric combination are, in fact, an exhaus- 
tive set of first order tensors. Using the properties of A and S listed above one can show 
that any other symmetric transverse projection can be represented as a linear combination 
of these structures and A^i,. Forming the trace-subtracted combinations we have, on writing 
out the tensors more explicitly. 



t 



-V 
2 



ApaAu/3 + A 



(2.5b) 



where r] is the shear viscosity. The parameter f/ is a P-violating transport coefficient referred 
to as the Hall viscosity. It is only allowed once parity is broken, and has been discussed 
previously for non-relativistic [16] and relativistic [7] fluids. If we denote small fluctuations 
of the spatial component of the velocity field by t;*, then unlike the normal shear viscosity 
which in flat space gives a response of the stress tensor T12 to (3i^;2+92'yi)5 the Hall viscosity 
gives a response of T12 to (^i-^i— 52^2)- 

There are four transverse vectors and four transverse pseudovectors which we can con- 
struct at first order in derivatives. These vectors and pseudovectors can be formed by pro- 
jecting VpT, V^(/i/r), F^^Uy = or u'^V^Up with either A^i, or T,p^. Since we have one 
transverse vector equation of motion, A^^V pT^^ = A^yF^fJp and one transverse pseudovec- 
tor equation of motion, TjpyS/ pT^^ = TjpyF^fJp, only two vectors and two pseudovectors are 
independent. We choose to drop the projections of u^V yUp. The constitutive relation for j'^ 
then takes the form, 



3^J. 



A 



aV^ + xeE^ 



XtV^T 



(2.5c) 



where Vp = Ep 



TApy\7'^(fi/T), and Ep is the electric field in the fiuid rest frame. 
We have written the constitutive relations in the Landau frame. However, a frame- 
invariant definition of the transport coefficients does exist and will be discussed in the next 
section. 
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scalars 


pseudoscalars 


transverse vectors 


tensors 




B — -If^'^^P?; F 
^ 2 "'/i-' vp 


Ul^ = F^"'u^ = EP 

= A/^-V,^ -^ = -}^ 





Table 1: Various independent first derivative quantities. The shear tensor a'^'^ was defined in ( 1.3c ) 
Pseudotensors and pseudovectors can be obtained from the vectors and tensors above through (3.5). 



3. Positivity of entropy production 

In the Landau frame the constitutive relations take the form given in (^]^) and ( |2.5D . In 
this section we wih study how the second law of thermodynamics restricts the coefficients in 
the constitutive relations leading to the results (|1.6|), (1.7), and (^) described in Section |^. 



At intermediate stages of the computation we will find frame invariant expressions for the 
constitutive relations. 

The assumption that the flux of entropy entering any compact spacelike region cannot 
exceed the amount of entropy produced in that region amounts to the existence of a current 
Js whose divergence is positive semi-definite, 

VpJ^^O, (3.1) 

with 

JM = ^Q^M + (j^^^ , (3.2) 

where sq is the entropy density given in ( [1.51 ). When there is no dissipation Js is conserved. 
The most general form of the entropy current in a 2-|-l-dimensional relativistic theory must 
take the form 

. . / all possible \ , , 

Js = "^icanon + ^'"S''^ '^'^''dicnt , (3.3) 



3- vectors 

where we will refer to 

^S^canon = ^0^^ " f T'^ - ^T>^^ , (3.4) 



with T^'^ and T'^ as in (|2.2| ), as the canonical entropy current. As we will see shortly, positivity 
of the divergence of the entropy current imposes non-trivial restrictions on both Jg and on 
t'"' and T/^ [1]. 

Our analysis closely follows [9,43]. We have described all possible first order transverse 
vectors, tensors and scalars in Section^ A list of independent transverse traceless symmetric 
tensors, transverse vectors, scalars and pseudoscalars is reproduced in Table |l| for convenience. 
Transverse pseudovectors U and pseudotensors S can be obtained from the above transverse 
vectors and tensors via 

iji^ = e^^'Pu^Up , §'''' = 1 (e'^^^tx^S/ + e'^'Pu^SpP) . (3.5) 
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transverse pseudovectors 




= e^^ ' UyV pT = —TUi — RoT 




— ^2 




= et^'Pu^Vpil. = IJP + ^ 




= ^eP^PF^p = + uPB 




= ep^py^up = -tj'^ + upn 



Table 2: An alternative basis for first order pseudovectors. In relating this basis to Table y we have 

used the equations of motion of ideal hydrodynamics to show that = RaTUt^, where 

Ra^ Po/{eQ + Po). 



We will refer to the first order terms which are independent under the equations of motion 
of ideal hydrodynamics as first derivative data. In this section we will find it convenient to 
use two alternative bases for the first order pseudovectors and pseudoscalars. The first basis 
is defined in ( |3.5| ) and Table ||. In addition to the basis of first-order pseudovectors spanned 
by BuP, QuP, and the Uj^^s we will also find it convenient to use a different basis of first-order 
pseudovectors, given in Table ^. In these bases the most general expression for the entropy 
current takes the form 

3 5 

= JL.non + M'^-u)nP + ^^Ut + E ^^^^ ' (3.6) 

1=1 i=l 

where the Vi's and the t'j's are (as yet) undetermined functions of fi and T. By including the 
pseudovectors V4 and V5 we have parametrized the longitudinal pseudovector contributions 
uPB and uP^l, respectively, to Jg. The entropy current (|3.6D contains all possible parity-even 
and parity-odd vector contributions both longitudinal and transverse to the velocity field. 
We note in passing that 



is a divergenceless vector for arbitrary a. If we add such a term to the entropy current ( p.6D , it 
will not contribute to entropy production but it will shift Di, and such that P5 — t- i>5 + d, 
z>i — 7- Pi — dxa and P3 — t- P3 — di±a. The two combinations of Pj's which are invariant under 
this shift are 8x1^5 + Pi and dRU^ + P3. Combined with P2 and P4, this means that there are 
only four combinations of the Pj's that a priori can participate in entropy production. 

The expression for the divergence of the entropy current in first order hydrodynamics can 
be written as a sum of products of first order data and a sum of genuine second order scalars, 

Vy Jl-l / products of \ _|_ / second ordciA f O Q\ 

^ S Vfirst order data/ \ scalar data J ' V*^*"/ 
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By second order scalar data we mean expressions which are second order in a gradient ex- 
pansion and cannot be decomposed into a product of first order terms. Equation ( |3.1| ) states 
that the divergence of the entropy current should be positive semi-definite for any flow which 
solves the equations of motion and with any background fields. Thus, all second order data in 
(|3.8|) should vanish and all first order data should arrange themselves into complete squares. 



Using the scalar equation u^V^T^^'^ + E^Ji, + fiV^J'^ = 0, together with ( p!^ and ([L5|), 
it follows that 



V.J°_ = - ( - ^ j - V, ) T^-' . (3.9) 



In other words, all the expressions in the divergence of the canonical part of the entropy 
current involve products of first order data. 

All second order data in the divergence of the entropy current must vanish. Therefore 
the coefficients, uq, Vi and Vi must be tuned so that no second order data appears in the 
divergence of the non-canonical part of Jg ■ Using (|3.9D it is not difficult to show that 



V, = + (1/2 - y) ^^^^^ + ^sA^^V^a,^ 



C, ^,CKV7 D I f products of \ 

U2U + (ji4 ^^^^^ j^,J 



where R^y is the Ricci tensor. Following an analysis similar to the one carried out in [9] one 
can show that all explicit expressions on the right hand side of (3.10) are genuine second 



order data. Since the genuine second order terms on the right hand side of ( |3.10| ) should 
vanish, we find that 

V2 = 1^0 = 1^1 = 1^2 = 1^3 = ^- (3-11) 

There are no additional constraints that arise from demanding positivity on a curved back- 
ground, so we will work in flat space from now on. 

We now require that the remaining data which contributes to the divergence of the 
entropy current appears quadratically. We allow the undetermined variables Vi , i = 1, 3, 4, 5 
to depend on fi = fi/T and T. To first order in the derivative expansion, 



UaJg — <-'a<^s canon / ^ 



4=1,3,4,5 



^{dT.Vi) + ^{dfl-V,) + i},{d-V,) 



(3.12) 
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Evaluating the right hand side of ( 3.12 ) we find 



a 

S canon 



-n{d-u) 

-B{d-u) 



T I I [OTVb + Vl) + —[ I (^M^5 + i^s) 



(3.13) 



+ U2-U3 [RoT{dTi>3 - d^Di) - dfiDi + RoT^dTi>A] 



+ Ui-U2 

where we have used the identities 



+ dpUi - dTi>3 - Tdr^A 



(3.14) 



CO 



derived in [9]. There is a redundancy in ( |3.13 ): Pi and 9^ only appear in the combinations 
9ri^5 + i^i and dfiu^ + i>3. 

We now evaluate the divergence of Jf,,„o„. Using dp) and we find 



"^S canon 



-A rf^" 



dPo 
den 



po 



dPr 



{Rou^.T^^'^ + T^) A.^U^ 



(3.15) 



2T 



In this expression, d-u is a scalar, A^^Un is a transverse vector, and is a transverse 



traceless symmetric tensor. Therefore, in (3.15), the first expression in the parentheses can 
be expanded in the basis of scalars, the second expression in the parentheses can be expanded 
in the basis of transverse vectors, and t^'^ can be expanded in the basis of transverse traceless 
symmetric tensors listed in Table |^: 



1 



den 



dpn 



f3{d-u) - (3nn - PbB , (3.16a) 



A^.A^.T^'' - ^A^^A^^r^^ 



al3 



ia 



(3.16b) 



and 



A°, {Rnu^T^''' + 



Tk3U^ - Tk3Ui . 



(3.16c) 
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The expressions on the left hand side of ( |3.16| ) are frame invariant. We refer the reader to 
[9] for further explanation and examples. The transport coefficients on the right hand side 
of (|3.16| ) are physical quantities which characterize the fluid and are independent of the frame. 
The relations between the parameters on the right hand side of ( 3.16| ) and the coefficients in 
the Landau frame are given in 

Requiring that the divergence of the entropy current be positive semi-definite leads to 
the following constraints: 

^Q=(|^) + {d;,i>5 + ^3), (3.17) 



in the scalar sector, 

Ki=0, K2 = 0, K3^0, K3GM, (3.18a) 

as well as 

Ki = T{dTi>4- Ro{dTh + , K2 = dpDi - RoidpD^ + i>3) , (3.18b) 
in the vector sector, and 

9^0, eeu, (3.19) 
in the tensor sector. A further condition is 

dph + h + T{dfiDi - drh) - T'^dTVi = , (3.20) 



which we used to simplify the expressions for ki and K2 in (|3.18| ). 

In the discussion below equation ( |3.7D we argued that only four combinations of vi par- 
ticipate in entropy production. Since we found that V2 = 0, three combinations remain: 
{Btv^ + vi), {dp_i>5 + Us), and z>4. Equation (|3.2Cl|) may be rewritten in terms of these three 
combinations, 

(dfih + h) + ^df, (t2 {drh + TdT {d-^v^ + i>3) - T^dTVi = . (3.21) 

We now parametrize the coefficients 1^4 and {djiv^ + P3) as 

TDi = MB. d-^D^ + D3 = ^d-^Mn-MB, (3.22) 
where Mb ^'iid ^-I'e arbitrary functions of ^ and T. Relation ( 3.21| ) then takes the form 

{T^drh + z>i) - TdrMn + 2Mn) = , (3.23) 
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which has the solution 



T'^idrv^ + 9i) = TdrMn - 2Mn + fn{T) , 



(3.24) 



where fniT) is undetermined. 

The Landau frame was defined by (|2.3| ), ( 2.5a ), ( p.5b ) and ( |2.5c| ). Inserting the Landau- 
frame expressions for the stress tensor and charge current into ( 3.16| ) we find 



1^1= XT, 
^1 = XT , 
Pb = XB, 



1^2 = XE. 
1^2 = XE- 

Pn = Xn- 



/3 = C, 

K3 = cr + RqTxt ■ 
K3 = 5- + RqTxt ■ 



(3.25a) 
(3.25b) 
(3.25c) 
(3.25d) 



Inserting ( |3.22| ) and ( p. 24]) into ( 3.17] ) and ( |3.18| ) and converting to the Landau frame ( p. 25 ) 
yields the expressions ( |1.8| ) for the x's as functions of AisiT, fi), Aifi{T,^) and fniT). 



4. Response functions 

The theory of linear response allows one to relate transport properties and thermodynamic 



we 



susceptibilities to limiting values of retarded correlation functions [51,52]. In Section 4T 
will discuss several properties of the retarded functions for fluids. These include positivity 
of spectral functions, covariance under time reversal, and connections to thermodynamic 



susceptibilities in appropriate limits. Then in Sections 4.2 and 4.3, we will show how these 



properties may be used to restrict the coefficients which appear in the constitutive relations. 
4.1 General properties 

Consider the expressions for the current and the energy-momentum tensor densities in the 
presence of external sources, 

J^{x) ^ {J^{x))a,, , T^'ix) ^ {T^''{x))a,, ■ (4.1) 

where x = (t,x), g = det{g) and g^^^ is the metric. The subscripts A and g indicate that 
these are one-point functions in the presence of a background gauge field and a metric 
g^iu = Vfiu + with ry^jy the mostly-plus fiat-space metric. When the background sources 
vary sufficiently slowly in space and time, the assumption of the hydrodynamic theory is that 



{J^^{x)) and {T^'^{x)) in (O) are given precisely by the constitutive relations ( |2.3D , (2.5), 
with the hydrodynamic variables eo, u^, and po satisfying the conservation equations (1.1). 
To obtain retarded two-point functions in flat space and without external flelds, we need to 
solve ( |1.1| ) to flrst order in A^ and /i^jy, insert those solutions into the constitutive relations, 
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and then differentiate (4.1) with respect to the sources. We define the retarded functions:^ 



GT{x) 



A=h=0 



A=h=0 



G'^^'^^x) 



6T^"'{x) 



A=h=0 



(4.2a) 
(4.2b) 



A=h=0 



Since we will be interested in the retarded functions in the equilibrium state which has B = 
and = 0, we need only solve the hydrodynamic equations of motion to first order in 
fluctuations in order to find and T^^^, to first order in background fields. 

Alternatively, the retarded functions may be defined in the canonical formalism, 

G^''^(x) = i0(t)Tr(^>[J^(x), rm), G>i^^^x) = ieit)TT{Q[T>^'^ix), J'^(O)]) , (4.3a) 



G''ji'"'{x) = ie{t)TT{e[r{x),T'^''m), G^"''^"(x) = i0(t)Tr(^,[r'^-(x), r^-(O)]) , (4.3b) 
where g is the grand canonical density operator specifying the equilibrium state of the system. 



The retarded functions defined by (4^) and (O) will differ by contact terms. See for example 
[55, 56] for a discussion. 

In the remainder of this subsection we will study general properties of the retarded func- 
tions defined in (|4.2| ). In the zero- frequency limit, the retarded functions may be determined 
by equilibrium thermodynamics. Indeed, the retarded functions at zero frequency coincide 
with the Euclidean time-ordered functions at zero frequency, and can be computed for exam- 
ple by the Euclidean functional integral method in the grand canonical ensemble. 

In a static equilibrium without external sources, the partition function Z in the grand 
canonical ensemble is given by 



Z[T, ^] = TV 



exp 



H uQ 
T T 



(4.4) 



where H is the Hamiltonian of the system and Q its total charge. In this coordinate frame 
the fluid velocity is n' = 0. When external sources are turned on, an interesting feature of 
the equilibrium partition function is that constant sources ^oo and hoi may be eliminated 
by a suitable redefinition of thermodynamic variables. To see this, consider the Euclidean 
theory. Let r be the Euclidean time with period /3 = 1/T. The chemical potential may be 
defined through the Wilson loop around the time circle 



/3 



drA-r 



(4.5) 



Suppose that the fluid is subjected to constant sources ^OO) and /loj. A constant contri- 
bution to ^0 shifts the chemical potential such that fi' /T' = fi/T + Aq/T. A constant /iqo 



^More formally, one can start with the generating functional in the Schwinger-Keldysh formalism [53, 54] 
and define the fully retarded functions by taking the appropriate variations as is done in, for example, [8]. 
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can be eliminated in a similar manner by rescaling r. Since the length of the time circle does 
not change under a rescaling, the periodicity of r must be shifted to /3' = — /ioo/2) to first 
order in /iqo- The temperature and chemical potential are then 



r' = r ( 1 + ^ 



/X =/i 



(4.6) 



to first order in external fields. Finally, a constant perturbation of the Minkowski metric /iQi 
may be removed by a coordinate transformation. The fluid velocity is invariant under this 
transformation, to first order in /iQj. Thus, the partition function in the presence of constant 
Aq, hoQ, and h^i takes the form 



Z[T, n] Aq, /loo, hQi] = Z 



T 1 + 



00 



M 1 + 



00 



+ ^o; 0, 0, 



(4.7) 



and the fiuid velocity is unchanged. Put differently, the thermodynamic behavior of a 
system with constant and small ^o, ^oo and hoi is equivalent to the thermodynamic behavior 
of the same system with zero Aq, Hqq and /loi but appropriately shifted temperature, chemical 
potential, and normalized three- velocity field. Retarded functions at zero frequency may then 
be evaluated in static equilibrium with n* = 0. By using the expression for the generating 
function on the right hand side of (|4.7| ) as well as the constitutive relations ( p. 2]) we find 



lim G°°'°(a;=0,k) 
k^o ^ ^ ' 



limG°°'°°(w=0,k) = r( , 



(4.8a) 
(4.8b) 



We will refer to the relations ( |4.8| ) as susceptibility conditions, and impose them as constraints 
that the hydrodynamic retarded functions evaluated later in this section must satisfy. 

The response functions defined in (4.3) are also constrained by their behavior under 
time-reversal T which is the basis for the Onsager relations [11, 12]. For two local Hermitian 
operators Oi and O2 which transform in a definite way under time reversal, @Oi@~^ = riiOi, 
the anti-unitarity of the time-reversal operator combined with translation invariance implies 

Gl{x) = ie{t)TT{Q[Oi{t,i^),Ojm) = i0(t)n,n,Tr (^>'[0,(t,-x), 0,(0)]) , (4.9) 

where q' = QqQ^^ is the time reversed density operator. We allow for the possibility that T 
may be broken in the microscopic theory by a set of real parameters such as fermion masses 
or the magnetic field, which we collectively denote as ha- The transformation T', which is T 
combined with ha — )• —ba, is then a symmetry, so that Q'{ha) = Q{—ba)- The Fourier transform 
of the retarded function must then satisfy 



G'^(w, k; ha) = niUjC^jliuj, -k; -ba) 
which will impose constraints on the transport coefficients. 



(4.10) 
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The canonical definition of the retarded function imphes 

ImG'^(w,k) ^ 0, (4.11) 

for oj ^ 0. When apphed to the Kubo formulas for transport coefficients, this condition can 
be used to argue that the viscosities t] and and the conductivity a cannot be negative. This 
agrees with the argument [1] that ry, and a contribute to entropy production and therefore 



cannot be negative. However, ( 4.11 ) does not constrain the sign of their parity-odd cousins 
fj, xn, and a. 

4.2 Parity-preserving hydrodynamics 

We first discuss the simpler case of fluids where parity is a symmetry of the microscopic 
theory. The constitutive relations in the Landau frame are given by (|2.3| ) and (^^) with the 
parity-odd terms omitted, 

T^^ = eouf'u'' + PoAf"" - CV/,n^ - r/^r^'^ , (4.12) 
= p^u^ + A^'^ [aVx + XeEx + XtVxT] , (4.13) 



where = F^"'uy, V^ = E^- TA^^Vifx/T) and 



Af_iaA,^IS + Aj^qA^^ — A^^Aal3 



V"n^. (4.14) 



While the constitutive relations are written in a particular frame, the correlation functions 
we obtain are frame invariant. They depend on the physical parameters defined in ( p. 161 ). 
In the parity-preserving case and in the Landau frame, these physical parameters are r], ^, a, 
XE, and XT- 

As we have shown in Section |3[ the positivity of entropy produciton implies that rj, 
and a are non-negative, while xt and xe both vanish. Alternatively, the same results follow 



from the properties of the retarded two-point functions described in Section |4.1| . We will see 
that the conditions 

C^O, o-^O, (4.15) 
follow from the positivity of the spectral function (4.11), while the conditions 



XT = 0, XE = 0, (4.16) 



follow from the susceptibility constraints (4.8). 



To compute the retarded functions we need to solve the linearized hydrodynamic equa- 
tions. The equilibrium state is described by a stationary and homogenous solution to the 
equations of motion ( |1.1D which has constant energy density eo, constant pressure Pq, the ve- 
locity field = (1,0,0), and vanishing sources. We solve the linearized equations of motion 
(1.1) in the presence of external sources h^i, and for 6T, du^ and 5u'^ which specify 



the linearized corrections to the chemical potential, temperature, and spatial components of 
the velocity field respectively. 
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The linearized hydrodynamic equations of motion may be written as 



DX 



(4.17) 



where X = [Sfj,, 6T, 6u^ , Su^) is a vector of hnearized hydrodynamic variables, D is a matrix 
of second-order differential operators, and 5 is a vector built out of the external sources hf^^, 
A^, and their derivatives. This system of differential equations may be Fourier transformed to 
obtain a set of algebraic equations for the Fourier components of the hydrodynamic variables. 
Taking the spatial momentum k in the direction gives 



D 



-lOJ 



V 



dfj. 

ikpo 




-ZOJ 



deo 



dT 

iksQ 




ikpo 

ik{eo+Po) 

- ioj{eo+Po) j 










(4.18) 



where sq is the entropy density, k = |k|, and all derivatives are evaluated at constant /i or 
T . The eigenvectors of the matrix D correspond to the hydrodynamic modes of the theory, 
whose dispersion relations may be obtained by solving for the roots of the determinant of D. 

Before presenting explicit expressions for the correlation functions, we pause to note 
that we are working with hydrodynamics to first order in derivatives, ignoring possible two- 



derivative terms in the constitutive relations.'* This means that equations (4.17) are valid 
only to second order in gradients of the hydrodynamic variables. In linearized hydrodynamics, 
possible two-derivative terms in the constitutive relations may be accounted for if we replace S 
with 5 + 0(5'^), and D with L' + C'((9^) where by 0(9^) we mean expressions whose combined 
powers of k and lo are greater than or equal to three. For example, the fluctuation of the 
transverse velocity bv? takes the form 



bu\S\ 



54 + 0(93) 



k'^T] - iw(eo+Po) + 0(&^) 



(4.19) 



We use parameterizations of this form to ensure that our results are free from these second- 
order corrections. 

To impose the susceptibility condition ( [1.8a| ) we compute the density-density function 



G^rf{uj,k) which can be obtained by solving the equations of motion (4.17) and using them 



to evaluate the one-point function Differentiating the resulting with respect to 



* Second-order hydrodynamics by itself is not a consistent efTective description of fluids, as the infrared 
effects of hydrodynamic fluctuations turn out to be more important than those due to two-derivative terms in 
the constitutive relations. See for example Ref. [57] for a discussion. We will ignore these fluctuation effects 
in this paper, which may be consistently done in the large- limit [58]. 
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Aq according to the definition (|4.2a|) , we find 



[CO, K) — — 



rdet(L>) 
dpo 



k\eo+Po) 



+ XrTpo 



where 



dpo dsQ dpo dso 



(4.20) 



(4.21) 



' ' dp dT dT dp 

is the determinant of the thermodynamic susceptibility matrix. One can expHcitly check that 
the 0{d^) terms do not contribute to the zero frequency, small momentum, limit of G^^{ljJ, k) 
and therefore we can reliably use ( ^30| ) to compute 



/ „ dpo dsQ 

r rO^O( n l-^ j. ^ ['^W ' P^W 
hm Gp [uj=0, k) = — h — ^ , — 



(4.22) 



The susceptibility constraint ( 4.8a ) then implies xe = . Similarly, a straightforward com- 
putation gives 



lim GD'^'^(a;=0, A;) 



(eo + Po)a 



'R 



{eo + Po)(T + TpoXt' 



dpo 
dT 



+ 



(eo + Po)TxT (dpo 
(eo + Po)(T + TpoXT \dp J J, 



(4.23) 



Using the second equation in ( [4.8a| ), one finds xt = 0. 



For the retarded function G 



12,12 



G 



R 

12,12, 
R ' 



one finds 



u, k=0) = -P + ir]UJ + ©(w^ 



(4.24) 



The positivity condition (4.11) applied to the spectral function of T^^ then implies that rj ^ 0. 
An identical analysis with the spectral functions of T^^ and shows that C ^ ^-nd a ^ 0. 
We have not found any other restrictions which follow from our consistency conditions. 

4.3 Parity-violating hydrodynamics 



We now move on to parity-violating hydrodynamics in 2-1-1 dimensions. Equations (2.3) 
and (|2.5| ) give the most general constitutive relations in the Landau frame in the absence of 



D 



parity as a symmetry. The linearized equations take the form ( 4.17 ) with 

fk^a-iu;^-k^^a + XT)-iuj^ ikpo \ 

ik{eo+Po) 

kHv+C) - iw(eo+Po) k^ (xn+v) 

-k'^fj k'^T] - iuj{eo+Po) ) 



-lUJ 



dtp 
dfj. 

ikpo 




-lUJ 



dtp 
dT 



iksQ 




(4.25) 
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Imposing the susceptibility constraints ( [4.8aD implies, via a computation identical to ( 4.22 
and ( 4.23| ), that xe = and xt = 0, exactly as in the parity-preserving hydrodynamics. 
To obtain Kubo formulas for the x's, it is useful to consider the combinations 



dpo 



dPo 
den 



7-00, C = T- RoT^\ 



(4.26) 



PO 



where Rq = po/(£o+^o)- A direct computation shows that at zero frequency, the 0{k) parts of 
C^, C* and are reliably evaluated in first-order hydrodynamics. We then find the following 
relations for the retarded functions at a;=0: 



liml(CV)/l(0,A;) 
fc-s>o tk 

lim^ {V'j^)R{0,k) 
k~^0 ik 



XB , 
-XE ■ 

0, 



liml(C°Oi?(0,fe)=XQ, 

fc-^o ik 

]im^{C^T^^)R{0,k) = -Tx7 
fc-s-o ik 

iimi(r*r°2)fi(o,fc) = o. 

fc-s-o k 



(4.27a) 
(4.27b) 
(4.27c) 



Note that in (4.27c) no summation over i is implied. For general k equations (4.27a) and 
( [4.27b| ) become the Kubo formulas ( |1.1C1D written in Section |l[ The last line will be relevant 
later in Section ^ 

The x's are also related by an Onsager relation. Using (4.10) to relate {C^C'^) to {C'^C^) 
and using the Kubo formulas (4.27) we find 

'dPo 



XB{ba) - RoXniba) 



XE{-ba) 



deo 



TxTi-ba). 



(4.28) 



PO 



If we assume that xe and xt are T'-odd (or alternatively, that xb and xn are T'-odd), then 
( [4.28| ) becomes 



XB - RoXn 



dP 



XE + 



dp 



TxT 



(4.29) 



where we have used the fact that the pressure is even under the symmetry T'. The same 
relation follows from the expressions (O) that we obtained for the x's by demanding the 
existence of a positive-divergence entropy current (there is also a differential relation between 
the x's which follows from (Il.8| )). 

We have not found any susceptibility or covariance conditions that relate or to other 
transport coefficients or thermodynamic functions. However, we do find Kubo formulas for 
them, 



lim 



lim 



lijj 



^ + XE: 



(4.30) 



which can be combined with limits of other retarded functions to give the Kubo formu- 
las ( 1.9a| ) and ( p..9b| ) presented in Section ^. Finally, the finite- frequency, zero momentum 
retarded functions ' , ' , and are the same as in the parity-preserving case. The 
constraints from positivity of the spectral functions are then unchanged, giving ( [4.15 ). 
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5. Thermodynamics 



The coefficients xe, XO) Xb and xt niay be evaluated from the zero- frequency Hmits of the 



retarded functions of the energy-momentum tensor and the current by (4.27). The retarded 
functions at zero frequency coincide with the Euchdean time-ordered functions at zero fre- 
quency, and the latter can be computed using the Euclidean functional integral method in the 
grand canonical ensemble. Therefore, it seems natural that the coefficients xe, Xn, Xb and 
XT should have a thermodynamic interpretation in terms of the derivatives of the partition 
function with respect to the parameters which leave the system in equilibrium. In this Section 
we will give such a thermodynamic interpretation to the functions Mb{T,^) and A^q(T, ^u) 
which determine the off-equilibrium entropy current, by treating B and 17 as parameters 
which characterize the equilibrium state. Equations (|1.8D then determine the coefficients xe-, 
XrJ) Xb and xt in terms of the derivatives of the pressure with respect to B and 17. 

Equilibrium states exist with non-zero magnetic field B, in which the pressure will depend 
on T, /i, and B. In addition, on compact manifolds one may have equilibrium configurations 
with non-zero vorticity Q.. Consider a uniformly rotating fluid on a flat disk with radius R and 
constant angular velocity satisfying uR <^ 1. Working to linear order in ujR the velocity 
field is given by 

= (l,-a;x^ ux^) , (5.1) 

where and are the two spatial coordinates, and the coordinate frame is assumed to be 
inertial, i.e. g^^ = = diag(— 1, 1, 1). The vorticity for this configuration is constant,^ 

-e'^'Pu^d^Up = 2uj . (5.2) 

The rotating fluid ( |5.lD satisfies the conservation equations (^]5) with constant energy and 
charge densities and exhibits vanishing V^u^, V^, and a^^ . To first order in derivatives the 
amount of entropy produced is determined by ( ^.13|) and ( 3.15| ), 



Since d^Js = for the flow described above we expect it to correspond to an equilibrium 
state. Other examples of equilibrium configurations with non-zero vorticity can be found, for 
example, in [59,60]. 

If vorticity is a parameter which characterizes an equilibrium state, the thermodynamic 
pressure will depend on ^, T, B, and 11, so that 

dP dP 

dP = sdT + pdfi + —dB + —dn, (5.3) 
oB as I 

e + p = sT + pn. (5.4) 



^ Instead of considering a rotating fluid on a disk, one could consider a rotating spacetime in wliicli tiie fluid 
velocity is u'^ = (1,0,0) which is related to ( ^.l| ) via a coordinate transformation. The value of the vorticity 
will, of course, be the same for both configurations. 
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To first order in B and ^l, the energy- momentum tensor and charge current in such an 
equilibrium state must take the general form 



T^"" = {e-eBB-enn)u^'u'' + {P-S:BB-xnn)A^'\ ^ ^ 

(5.5) 

= {p-rsB -rnn)u^' . 



where cb, en, xb, xq, tb, are functions of T and /i. Since ( |5.5D describes an equilibrium 
configuration, the quantities cb, en, xb, Xfi, tb, ro are all measurable. In the remainder of 
this section we obtain the most general expression for these parameters within the framework 
described in the previous sections. Indeed, the expressions (5.5) fall into the class of energy- 



momentum tensors and currents described by ( |2.2| ), where we now view the terms in r^'' 
and T'^ proportional to B and Q as equilibrium quantities. Therefore, we are free to use our 
results from Sections ^ and ^ specialized to equilibrium configurations with non-zero B and O 
in our analysis. 

Since we are assuming that B and 17 label equilibrium states, the Euclidean partition 
function Z will depend on fi, T, B, and $7. Recalling footnote ^, it is convenient to work in a 
coordinate system where = in equilibrium. To linear order in fluctuations, the vorticity 
is $7 = diu^ — d2U^ + e^^dihoj. However, for this equilibrium state and in this coordinate 
system, diU^ is higher order in k and the vorticity reduces to O = e^Wihoj. In this case, since 
small shifts of the magnetic field B = e^^diAj and the metric perturbations e^^dihoj do not 
affect the radius of the Euclidean time circle or the Wilson line around it, T and /i should 
not change to first order in B and 0. Thus, using the same formalism as the one described 



in Section which led to (|4.8| ), we expect that 



lim ^GfiO, k) = ^-rB, hm ^Cf^O, k) = ^ - cb , (5.6a) 

k^o ik oB k^o ik oB 



where r^, r^, e^, and were defined in (p?^). 

Now we recall the Onsager relation ( [4. 281 ) obtained at the end of Section |^. Combined 
with the expressions for the x's in ( |1.8D , obtained via the entropy current, we find that the 
combination {dPo/dpo)xE + idPo/d€o)TxT is odd under the symmetry T', which is time- 
reversal T combined with sending the T-violating parameters ha ^ —ha- In what follows we 
assume that both xe and xt are T'-odd. 

Applying the Onsager relations ( 4.10| ) to the Kubo formulas ( 4.27b ) for xe and xt we 



find 

XE = lim 1 (g°/(0, k) - i?oG?^°'(0, k)) , TxT = Hm 1 (cf^O, k) - RoGf'\o, k) 

k—i-O ik \ J k-^O ik \ 

(5.7a) 
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The Kubo formulas for xb and xn in ( 4.27a ) are also linear combinations of the same retarded 
functions, 

1 fdPono, 



XB = lim 



fc-i-o ik \ dpo 



-Gj^ {0,k) + —Gj^' {0,k) j , 



1 fdPo^O,02rr. ,x , ^^0, 



Xn = lim - ( (0, k) + ^Gf'\o, k) ) . 

fc-i-o ik \ opo oeo 

-,0,2 ^00,2 ^0,02 „^ , ^00,02 



(5.7b) 



Solving (^j) for , G^ ' , , and G^ ' and using (^^) we obtain 



fc-^oifc ^ ^ ^ 5/i 



lim^ TtG^^ (0, k) = , {51 



Comparing ( |5.8| ) to { p.6\ ) we find 

■Mb 

OMb _ dp_ 
dp ~ dB 



BMb ^ OMb de 

+ p^ -^^ = dB-'^^ (^-^^^ 

- TB , (5.9b) 



By defining = vq — M.b and = — A^q + fn, equations { p.9\ ) can be brought to a 
more symmetric form. We can solve for Mb, -Mn, gb and en in terms of the functions vb 
and rn, 



eB = T^-^ (^i p rBdp^ + ^rs - T^^T) , 
en = y fod/i^ + pfu - T'^q'^iT) , 



dP 



(5.10) 



^B=Q^-J rBdp + TqB{T), 

dP ff" 
A^n = ^ + / rndp + TqniT) , 

where qB{T) and qn{T) are undetermined functions. 

The coefficients xb and can be determined from the retarded functions 



lim ^Gt^''(0, k) = ^-XB, lim ^G'^^\o, k) = ^ 
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with no summation over i. Using ( 4.27c| ) and ( ^.llD we find that 

dP dP 

If B and 17 are indeed equilibrium parameters, we should find that in an equilibrium 
state characterized by T, /x, B and 0, 



(5.12) 



-So 



ds 



B 



Going back to our expression for the entropy current in 
T'^ and r'^^ from (2.2) and (p.5|), we find that (5.13) implies 

dP 



TD4 
TPs 



dB 
dP_ 



+ {en - urn) . 



(5.13) 

and using the expressions for 

(5.14) 
(5.15) 



Remembering the definition A4b = Tu^ and comparing 

( CT) with (|5l0| ) we obtain 



driqBT) 



(5.16) 



from which we conclude that tb = fBifi) and qB{T) = Qo/T, with constant go- Defining 
hsilj) = Qo — ^B{tJ')dfi we find, for $7 = 0, that the equilibrium energy-momentum tensor 
and current take the form, 

dP 



rpflV 
J" 



(e + {fih'B - hB)B) 



P 



dB 



B A^" 



{p + h'^B) 



and 



M 



B 



dP 
dB 



+ hBin) 



(5.17) 
(5.18) 

(5.19) 



The result ( 5.17 ) with Kb = is the canonical expression for the stress tensor and current 
of a fluid at non-zero B and 17 = 0, see for example [4,61]. Following the literature, we set 
Hb = 0, although we have not found a constraint that eliminates it. Note that Hb = implies 
that cb = rB = 0- Then the only effect of the magnetic field is to shift the equilibrium value 
of ^Afj^iyT'^'^ away from the thermodynamic pressure by the term —{dP/dB)B. We will refer 
to this difference as a "subtraction". 

Unfortunately, a similar analysis involving the entropy current in equilibrium cannot be 
carried out for the coefficients rn and en- As emphasized below equation (3.7), the definition 
of the entropy current is inherently ambiguous. As a result, fixing through ( |5.15| ) without 
additional knowledge of ^3 or Pi does not provide additional constraints on the thermodynamic 
response parameters. However, based on the similarity of the equations for fn, en and Mn 
in (^) with those for r^, and A^b, we conjecture that 

dP 

rn = en = 0, Mn = ■ (5.20) 
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(5.21) 



This leads us to the equilibrium constitutive relations ^ 

T^- = (e - {Mn - fn)^) u^'u'' + [p - (MbB + Mn^)) A'^" , 
J^' = {p- Mb^) u^' , 

together with 

dP dP 



These relations determine the 'magnetovortical' frame ( 1.13 ), and we will discuss the above 
results in Section |^ 



6. A holographic model 

In this section we obtain explicit expressions for the constitutive relations, transport coef- 
ficients, and thermodynamic response parameters of parity- violating hydrodynamics in 2+1 
dimensions in a relatively simple model. This is made tractable by the AdS/CFT correspon- 
dence. Our holographic setup is described by the action 



(6.1) 

where 6[ip\ = 6q + 6iip + 0{if'^) is a function specifying the parity-violating axion coupling, 

the potential satisfies ^[0] = V'[0] = 0, and 5bdy is an appropriate boundary contribution 

to the action. We will denote the coordinates as = {x^,r), where r foliates the spacetime 

such that r — )• oo is an asymptotically AdS boundary, and are coordinates on a constant 

r hypersurface. We denote the completely antisymmetric Levi-Civita symbol by e"'^^'^ where 
g0123 _ ^ 

The AdS/CFT correspondence [62-64] putatively relates the gravitational theory defined 
by (|6.1[) to a strongly coupled 2 + 1-dimensional theory which can be thought of as living on 
the asymptotically AdS boundary of the spacetime located at r — )• oo. Equation ( |6?l| ) defines 
a simple two-derivative action which can describe parity violation in the charged sector of the 
dual field theory. There are a number of string-theoretic embeddings of ( |6.1| ), including the 
worldvolume theory of probe branes dual to fiavor multiplets [65] , as well as truncations of 1 1- 
dimensional supergravity [66]. A similar holographic model, with parity violation introduced 
through a gravitational axion rather than an electromagnetic axion, has recently been studied 
in [10]. 



®Using ( 3.16a| ) and (3.21:), we find that the parameters xb ~ dP/dB — Mb and xq — dP/dO. — Mn 

are related to xs and xn in the Landau frame (L2) by is = Xs + (|§ - H^fI ~ l^ls)' in = + 

iLQi ^ A/( f ^ _L iLQi A/t ^ ( sp aPo at app dp \ 
(Mn - fn) + -b^Mb + [gn ' al^ an ' a^ an ) ■ 
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The action (6J) leads to the equations of motion: 



R, 



R 



ab 



where we defined 



Tab 



da 
da 

1 

2 



-Qab 



ab 



L2 



9ab + Tab : 



-gF 



J)cde 



87r2 



drOF, 



de ■ 



(6.2) 



Fab Fed ) 



1 



gabF"^ + da^pdb^p - -9ab{dipf - gabV[if\ 



(6.3) 



which is not to be confused with t^'^ from (2^). In the remainder of this Section we choose 
units in which L = 1. 

The AdS/CFT correspondence relates solutions of (|6.2|) to expectation values of operators 
in the dual 2 + 1-dimensional field theory. In particular, the gauge field is dual to a 
conserved current J'^ whose expectation value is given by 



{J^{x)) 



1 



5S 



1 



lim 



2k^ 



+ 



r) 



(6.4) 



where gB is the background metric of the dual field theory. Similarly, the scalar field (f is dual 
to a neutral scalar operator O^p. The conformal dimension A of O^, satisfying 3/2 < A < 3,^ 
is related to the mass of through V"[0] = A(A — 3). We will break conformal invariance 
by turning on a source term = A'^~^ dual to O^, with A characterizing the strength of 
the source. In other words, we will deform the field theory action by a j_a J d^XyJ—ggJ^pOip 

2 

term, where the prefactor (A — |)^^ was introduced in [67] in order to have a smooth A — )• 3/2 
limit. Near the boundary, we find according to [67] that 



A3 



-A 



l + 0(r-i)) + 



3 
2 



^ (l + 0(r-i)). 



(6.5) 



A-l ^ — ^ A 

We will restrict our attention to the high-temperature regime T ^ A, T ^ /i which evades 
potential low-temperature instabilities, and work perturbatively in small /i, J,^, and Vt. The 
expectation value of the energy-momentum tensor of the boundary theory can be computed 
along similar lines to [68,69]. However, in this work we will only need the trace of the 
boundary theory which can be computed from using 



A 



a3-^(o, 



(6.6) 



The Ward identity ( |6.6| ) can be derived holographically as in [69], or by a field theory com- 
putation as in, for example, [70]. For particular values of A there may be an anomalous 
contribution to the trace of the stress tensor associated with a matter conformal anomaly. 
We study values of A for which this anomaly is absent. 

^Unitarity implies that A > 1/2, while is a relevant operator if A < 3. We focus here on A G (3/2, 3) 
for computational convenience. We expect that our results may be analytically continued to 1/2 < A < 3. 
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6.1 Equilibrium configurations 

In this subsection we construct equilibrium solutions to ( |6.2| ) to first order in the magnetic 
field and vorticity. Various solutions to (6.2) with ip = and non-vanishing magnetic field 
and vorticity do appear in the literature [71-74]. However, a more useful starting point for 
our analysis is the solution described, for example, in [75] (see also [76,77]), 

3 

gabdx'^dx'' = 2dtdr + r\-f{r)de + dx^) , /(r) = 1-"^ + 0{f,\ JI) , 



Aadx'' =[^Ji-^)dt + 0{4), (6.7) 



where P\{z) is a Legendre function of the first kind. The temperature of the boundary theory 
is given by 

T=^-^ + 0{^,^Jl), (6.8) 



the chemical potential is = q/rn, and the charge density is p = qjlK^ . 

We would like to turn on a small magnetic field and vorticity such that the boundary 
theory stays in equilibrium. Consider the ansatz. 



dgabdz^dx^ = 2Qr^f{r)x^dtdx'^ - Qx^drdx'^ + 0{Ji 



^)x^dx^^O(J^.). ^'-'^ 



SAadx'^ =[B + x^dx^ + 0{Jl) , 



which describes a configuration in the boundary theory with non-zero vorticity and non- 
zero magnetic field B. When 9[ip] = 0, ( |6.9D solves the equations of motion (6.2) to linear 
order in B and Cl. We may now solve ( |6.2| ) for 9[ip] ^ by further perturbing the solution in 
(|6.9D to linear order in B and il. We find 



47r2 

Sip = — 



+o(^^ 4), 



(6.10) 

K'^Oifirn r dr' r' dr" ^ir") ( ^ firHn\ ,^,0^ . 



where 61 is defined below ( |6.lD , and oq is a constant, chosen so that 5Aq vanishes at the 
horizon, 

""^ dr r , , „ f „ qn' 



T . 



"0 = / -2 / dr'^'ir') i3 + ^ . (6.11 



In obtaining ( |6.10| ), we have applied equilibrium boundary conditions to the equations of 
motion. Namely, regularity at the horizon and normalizability at the AdS boundary. This 
last condition corresponds to holding the sources of the boundary theory, fi = lim^^oo ^0 and 
= A^"'^, fixed. The perturbations SAq and dip also source metric corrections, whose effect 
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on susceptibilities is suppressed at order /j.'^J^. The magnetic field and vorticity appear in 
the combination B + q^/r, which is a consequence of Q, appearing in (|6.9|). 

Using ( |6.4D and ( |6.5D and integrating by parts we find that the boundary theory current 
and scalar expectation values are shifted by 



+ ^ + 0(M=,(.2B,,.2n,J^), (6.12a) 

*(0,) =A-' (a - ^ £ ^ (b + fl!^) + 0(,^B. . (6,12b) 
In what follows it will be convenient to define 

where ^{r) is given by (|6.7|). Using ( |6.6D we infer that, at B = = 0, 

^ = ^^ + ^2+ 4) ' ^ = 2^-^^3 - ^^2 + J^) , (6.14a) 

= (3 - A)m<I>2 + J^) , = (3 - A)^rH<l>3 + 0{^\ 3%) . (6.14b) 

We are now in a position to compute all the susceptibilities. In particular, using the expres- 
sions J° = p - TWfifi, TOO = e + {fn{T) - Mn)n, T^^ = T^^ = P - MbB - MqQ and 
dP = sdT + fidp + MndB + Mud^l from Section |, we find 



ll = ^2 + 1^ + 0(/.2, J2) , ll = - 3)/^'&2 + 0(/x2, J2) . (6.15a) 

Integrating the Maxwell relation dp/dB = dMB/dp leads to 



^=MB=P<^2 + f^ + cb{T) + 0{p\ 4) , (6.15b) 

where cb is the magnetization density at p = 0. However, cb vanishes due to the symmetry of 
the model under charge conjugation C. Having obtained A^_b, we find the vortical variation 
of p, via = p — M. B^, to be 



^ = 2^,rH^^ + ^ + 0{^?, 4) . (6.15c) 



We then integrate the Maxwell relation dp/dO, = dA4n/dp to find 
OP lP'O 

— =Mn = P^rH<^3 + + cn(r) + Oip^ J^) , (6.15d) 

where cq is dP/dO, at ^ = 0. We now point out that at S = and p = the perturbation 
is an exact solution to ( |6.2D to linear order in It follows that the on-shell bulk action 
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density, equal to the pressure of the boundary theory, is unchanged at this order so that 
dP/dQ vanishes at /x = 0. This ehminates cq. Using the expressions for T^^ and T**, we find 
the final susceptibility to be 

£ = (A - 2)^V^ci>3 _ /^(T) + Oif,^ 4) . (6.15e) 

We can determine fn{T) by using 

TdT^ = rndr^^ = {A - 3)ip - rif' , (6.16) 
which follows from (|6.7]), to show that the magnetic and vortical variations of e are 

^ = T^ + ^^-A1^, (6.17a) 
de ^dMn dMn ^rtvk^ 



at the order to which we are working. The thermodynamic relations ( |6.17a ) and ( |6.17b ) 



fohow from ( |lTT|) , dP = sdT + pd^i + MsdB + MndO-, and (p^) , e + P = sT + up, only 
when fn{T) = 0, implying that /q vanishes. We emphasize that in obtaining ( 6.15| ) we 



used ( 5.21D without deriving it explicitly in the holographic model. However, we note that 



( |6.17| ) represents a nontrivial consistency check of ( |5.21| ). 
6.2 Fluid-gravity correspondence 

The fluid-gravity correspondence [33] provides a method to explicitly compute the hydro- 
dynamic constitutive relations of the boundary theory. In this subsection we will use it to 
check the consistency of the parity- violating constitutive relations described in ( |1.13| ), ( 1.15 ) 



and ( |6.15| ) . We begin with the background (|6.7D boosted to a velocity u'^ normalized so that 
u^rj^yu'^ = — 1, with rj the Minkowski metric. In what follows we define n^j = •q^yu'^ and 



gabdx^dx^ = -2u^dx^'dr + r\-f{r)u^Uy + A^^)dx^dx'' + 0{p\ J^) , 



Aadx'' = ^u^dx'' + 0{p', J^) 



(6.18) 



with /(r) as in ( |6.7D . We would like to retain the distinction between a chemical potential 
and an external background field in an out-of-equilibrium configuration. For this reason, the 
gauge field in ( |6.18 ) has been shifted relative to the one in (pTz]). When n^, T, and p are 



constant, the profiles (6.18) solve the equations of motion 



Consider an extension of the solution ( |6.18| ) where p, T and vary slowly as a function 
of the x^ coordinates, and, in addition, an external electromagnetic field whose magnitude is 
of the order of the gradients of p, T and u'^, is turned on. This extension of ( |6.18 ) no longer 
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solves the equations of motion (|6.2D , so the bulk fields need to be corrected. Since we have 
chosen u'^, T and /x to depend slowly on the spacetime coordinates, we can correct the bulk 
fields order by order in derivatives of n'^, T and fi. We write 



^ = ^(0)+^{i) + ... 



+ ^(1) + . . . 



(6.19) 



where g^^\ A^^\ and 99^^) refer to ( |6.1^ ) with slowly varying T, fi, and n^. See [33-35] for 
further details of the derivative expansion. 

We parametrize the corrections to the metric and gauge field by 



g'^^dx°'dx^ = -2gi{r)Uf,dx^dr + r'^{go{r)u^Uy + 7^^,, + 7,,n^ + ir^y + gi{r)Af,^)dx^'dx'' 



A'^^Ux" = -ao(r)n^dx^ + a^{r)dx^' + A^^^dxf" , 



(6.20) 



where Oif^{r) and 7^(r) are transverse to u^, TT^^i^) is symmetric, traceless, and transverse 
to u^, and j4jj^(x'^) is an arbitrary electromagnetic background gauge field, to first order in 
derivatives. To first order in the chemical potential, we can decouple the equations for 7^(r) 
and a^(r). We find 



4^2 



bg,p 



dpq + qu^'daUi: 



(0)' 



(0) 



+ 0(/z2,j2), (6.21 



with q = r/f/i and £'^g = F^^gU^- Equation ( |6.21 ) may immediately be integrated to obtain 
Op and, via (|6.4| ), the transverse component of the boundary current. We now translate this 



result into the constitutive relations given in ( p. 5c ). To this end, we need to relate u'^dp^Up 
to A'^'^duT, A^'^dufJ' and Ej^^ using the equations of motion of ideal hydrodynamics. Energy- 
momentum conservation of the boundary theory is encoded in the Einstein equations, as 
described in [78], yielding 

u o^u'^ = ^ , (6.22) 

which is just the transverse vector part of the ideal hydrodynamic equations. After using 
q = rnfJ-, (|6.8|) , ( |6.15D and ( 6.16| ), we find that the constitutive relations of the holographic 
fluid take the form ( |1.13| ) with 

dp 



dB 



— +o(/i^J^), 



(6.23a) 
(6.23b) 



where the susceptibilities are given in ( |6.15| ), and the Hall viscosity vanishes, fj = 0.^ The 
constitutive relations with (|6.23|) match our earlier predictions (1.8) once we recall that 



*The linearized equations of motion for the tensor modes of the metric /112 and (/in— ^22), which follow 
from (6.1), decouple from the equations for Aa and </p at fc = 0. The Kubo formula (4.3C) then implies that 
?7 = is an exact relation in this model. 
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Mb = Oifi), Mn = ©(^u^), dp/dn = 0{n), and de/dQ = 0{^i^), Rq = 0{fi), so that 
their contributions to xe and xt are ©(/u^) and 0{^^), respectively. We also note that at 
this order the off-diagonal conductivity is given by 

a + XE='^^ + Oi,-), (6.24) 

which is reminiscent of the membrane paradigm [79]. 

We may compute the trace of the stress tensor by solving for the correction to the scalar 
field and using (|6.6| ). At high temperatures the scalar equation of motion becomes 



(/rV(^) ')' - r^A(A - 3)^« = [b^, + ^) + (-rr ) , (6.25) 

where we have used the definitions Sbg = and = —e^'^fuydi,Up. The parity- 
violating part of the solution for is the same as that for Sip, (|6.10[) , as described in 
the previous subsection. Therefore, the shift in the trace of the energy-momentum tensor 
coincides with that described in ( |6.14b ), 



where the omitted terms are independent of B and and we have used the equilibrium 
susceptibilities ( |6.15| ). It is satisfying that the magnetic and vortical subtractions computed 
in equilibrium match the dynamical result as expected on general grounds. 



7. Discussion 



In this work, we have presented a general framework for first-order relativistic hydrodynam- 
ics in 2 + 1 dimensions, applicable to systems with parity violation. We have attempted 
to describe our computations and results in as simple a form as possible, avoiding unneces- 
sary interpretation, with the results summarized in Section |l[ In this section we include a 
discussion of a few properties of this hydrodynamic framework. 

• Thermodynamic and hydrodynamic response. The Kubo formulas computed 
in Section |^, and summarized in Section |l], naturally split into zero momentum (or 
hydrodynamic) correlators and zero frequency (or thermodynamic) correlators. The 
parameters r/, i^, a, fj, and a obtained from zero-momentum correlators are transport 
coefficients, while xe, Xt, Xb and xn obtained from zero-frequency correlators are 
thermodynamic response parameters. The difference between hydrodynamic and ther- 
modynamic response parameters has been discussed previously in the literature. See, 
for example, [8, 33, 80] for a discussion of thermodynamic response parameters in sec- 
ond order hydrodynamics in 3 -|- 1 dimensions, or [36] for examples of thermodynamic 
response parameters in first order parity-violating hydrodynamics in 3 -|- 1 dimensions. 
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Curiously, the Kubo formulas for the anomalous 3 + 1-dimensional transport coefficients 
and 4, defined through the Landau frame charge current 



4+1 = puf" + aVf' + ^33" + Cw'', (7.1) 

with B^^ = ^e^^^f^UyFpfj the magnetic field and w^^ = e^^^P'^UyV pUa the vorticity, are 
tantalizingly similar to the Kubo formula ( 1.1C| ) for XE and XT- 



Thermodynamic subtractions. In Section H we obtained the equilibrium expressions 
for the energy- momentum tensor and the current ( ^.21 ), 



T^- = (e _ i^Mn - fn) u'^W' + (P - (MbB + Mnn)) A^^ , 
JA* = {p-MB^)u^'. 

The above structures can be given a physical interpretation. First, the magnetic sub- 
traction to the thermodynamic pressure can be understood as being due to the force 
exerted by B on the boundary currents of the fluid elements [61] . Second, a rotating 
fluid with non-zero magnetization must also have non-zero polarization (see for exam- 
ple [81]), thereby inducing effective "bound charges." In [4] the latter effect manifests 
itself as an extra correction to the current which takes the form = pu^ + d^M^'^. 
Rewriting M^'^ = —Mb^'^'^^u^, taking p and T to be constant and using the equa- 



tions of motion (1.1) to show that the transverse part of e^'^^d^Up vanishes, we obtain 



= {p-Mb^)u'' as in (| 

Restrictions on the constitutive relations. Had we known the equilibrium stress 
tensor and current at non-zero B and O, we could have obtained the susceptibility 



conditions (|5.§| ) without requiring any input from the entropy current. Applying (5.8) 



and the Onsager relations ( 4.10 ) to the Kubo formulas for the x's ( 4.27 ) then determines 
■Mb, -Mq, Xe, XT, xb and xq to take on the same values we found by demanding 
positivity of entropy production. It is interesting that the constraints on hydrodynamics 
from correlation functions then precisely match those obtained from the entropy current. 
In each case the x's are functions of dP/dB and dP/dfl and the integration constant 
/n(T). Note that dP/dB can be interpreted as the magnetization density while, as 
discussed below, dP/dVl can be interpreted as half the angular momentum density of 
the equilibrium state. 

Equihbrium states with vorticity. The possibility of equilibrium configurations 
with non-zero vorticity O was discussed in Section ^. The example we gave involved 
a limit in which the vorticity was taken to be small, of the same scale as gradients 
of thermodynamic variables. In this limit, there is a solution to the hydrodynamic 
equations in which the fluid is rotating with constant angular frequency oj and has 
vorticity O = 2a;. We point out that a fluid rotating in a flat background is equivalent, 
via a diffeomorphism, to a non rotating fluid moving at constant velocity, in a rotating 
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background. If we denote Bq = e^Wihoj, also known as the gravitomagnetic field, we 
find using the aforementioned diffeomorphism that Bq = 2a;. This relation has been 
used recently in a study of response in topological superconductors and superfluids [82]. 
Equilibrium states in the presence of gravitomagnetic fields were considered previously 
in [15]. We observe that diffeomorphism invariance implies that equilibrium states 
characterized by either a constant angular frequency or a constant gravitomagnetic 
field are equivalent to states with non-zero vorticity (see also [36]). While the vorticity 
is diffeomorphism invariant, Bq and u are not. 

There also exist equilibrium configurations where does not need to be small. One 
example is given by a fluid on a sphere [59] where Q, = 2lo cos 0/(1 — sin^ 9), and 6 
is the polar angle on the two-sphere. Another candidate example is an extension of our 
solution from Section ^, where u'^ = 7(1, 0, co), 7 = 1/ y^l — (cjr)^, and we are working 
in a polar (t, r, 6) coordinate system. In this conflguration we find that Q = 2uj^'^ while 
the shear and divergence of the velocity field vanish, afj_^ = and V ^u^ = 0. See [60] 
for additional examples. 

Vorticity and angular momentum. In an equilibrium configuration where either 
the spatial fluid velocity or angular velocity are constant but non-zero, the partition 
function (14.41) takes the form 



Z[T,/i,n^] = Tr 



expi Jd^xiu^T'^^ + fij') 



(7.3) 



Inserting the flat space rotating solution ( p.lD from Section g which has non-zero vortic- 
ity ri, together with the definition L = f d?xeij x'^T^^ for the total angular momentum, 
leads to 

Z = Tr exp (-— + — + — 1 . (7.4) 
We observe that ( |7.4| ) describes an ensemble with non-zero total angular momentum. 

Ferromagnetism and Ferrovorticism. When Q characterizes equilibrium states, 
an interesting feature of P-violating systems is that the response of the pressure to 
vorticity, Ain, and to a magnetic field, A^_b, can be non-zero when Q=0 and B=0. 
While a non-zero value of magnetization density Mb = dP/dB at i?=0 is associated 
with ferromagnetism, one may term a similar phenomenon of non-zero A^n = dP/dO. 
at 17=0 "ferrovorticism". For the rigid rotation states (5.1) discussed in Section ^ with 
constant vorticity 17 = 2a;, the form of the partition function (^^) shows that we have 



Mi 



dP 



1 dP 

2 



1 



(7.5) 



aj=0 



where {i) = {L) /V denotes the density of the total angular momentum at zero rotation 
frequency. In parity- invariant fluids, the total angular momentum in the parity-invariant 
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thermal equilibrium state with a;=0 would have to vanish. However, when parity is not 
a symmetry of the microscopic theory, there is no reason why {L) should vanish at a;=0 
(just as there is no reason that the total charge {Q) should vanish at /i=0 when charge 
conjugation is not a symmetry). The holographic model studied in Section ^ indeed has 
a non-zero value of dP/dQ at Q=0, given by (6.15). 



Recently, Ref. [7] argued for a relation between the angular momentum density and the 
Hall viscosity at T = 0. It is worth pointing out that our holographic model of Section ^ 
has a non-zero value of Ai^ = dP/dO, at T 7^ 0, while the Hall viscosity fj vanishes 
identically. 

The ofT-diagonal conductivity. The total off-diagonal conductivity axy, defined 
through = axyE^ , is given by 

(Txy = ^ + XE , (7.6) 

where xe is a thermodynamic response parameter and a is a transport coefficient. In the 
holographic example presented in Section ^, a 7^ 0. We anticipate that a will contribute 
to an anomalous Hall effect [14] in parity-violating systems, and to an off-diagonal heat 
current, = r°* - fiJ\ at non-zero chemical potential. We point out that the total 



conductivity in (7.6) contains contributions from both transport and bound currents. 
We plan to study these contributions in future work [83] . 

The fluid-gravity correspondence allows an interesting interpretation of the off-diagonal 
conductivity in terms of properties of the black hole horizon. The fluid-gravity analysis 
of Section ^ can be extended to include a dilaton, e{x'^), coupled to the gauge field 
via a e~'^{x)F^yF^^^ term which replaces the canonical kinetic term. Such an extension 
will lead, at zero chemical potential, to the DC longitudinal conductivity axx = o = 
Xje^irH) along with the off-diagonal conductivity Oxy = 5" -|- xe = ^('^_ff)/(87r^).^ Since 
only the horizon value r = rn oi the axion and dilaton enter the expressions for the 
conductivities this result may be suggestive of the membrane paradigm. In spite of 
this possible interpretation, one can show that at non-zero chemical potential the off- 
diagonal conductivity cannot be written in terms of horizon quantities. 
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